1. IntroductionDiamond (C) has attracted wide interest as an extreme semiconductor material for next-generation power and photoelectronic devices due to the wide bandgap, the highest thermal conductivity, highest breakdown field, highest carrier mobility, and highest carrier saturation velocity in all the available semiconductors.[1] The practical application of semiconductor diamond electronic devices relies on the reduction of the electrical resistivity of diamond by doping impurities and controlling the interface formed between diamond and other material. Despite the lack of shallow dopant at this moment, the hydrogenated diamond surface exhibits a unique p-type surface conductivity with a hole density of around 1013 cm−2.[2] Such a surface conductivity has opened the avenue to the development of metal-semiconductor or metal-oxide-semiconductor field-effect transistors (FETs), which showed exciting performance with either high cut-off frequency or capability operating at elevated temperatures.[3,4] Another important structure is the hetero pn-junction, which is the basis of the development of deep-ultraviolet photoelectronic and power devices.[5] Among all these heterostructures, the interface properties at the junctions of diamond with other materials determine the ultimate device performances.
One promising junction is the AlN/C heterostructure, in which both AlN and diamond are wide-bandgap semiconductors. The AlN/C heterostructure is expected to lead to a variety of novel semiconductor devices such as deep-ultraviolet (DUV) light emitting diode, DUV photodetector, and high-power FET. Experimentally, the implementations of AlN/C pn-junction diode and field-effect transistors have already been demonstrated.[6–8] Although single-crystal AlN epitaxial layer on diamond could be achieved, the existence of a large number of interface states is anticipated due to the large lattice mismatch between the hexagonal AlN and the cubic diamond. The high concentration defects degrade the device performances such as the occurrences of sub-bandgap emission and the low efficiency of band-to-band recombination in pn-junction diode and high leakage current in the electronic device. Therefore, the understanding of the interface characteristics of the AlN/C heterojunction is in demand for the development of the next-generation optoelectronic devices operating at DUV region and electronic devices for high-power and high-frequency applications.
However, few experimental values[7] or theoretical results[9,10] are available in the interface characteristic parameters of the AlN/C heterostructure. The reported calculations are only limited to the band offset of polar plane AlN(0001)/C, and the investigations on the electronic properties at the interface of semi-polar plane AlN(
)/C heterostructure is still lacking. In order to fully realize the electronic properties such as the electronic structure, band-offset, and charge transfer at the AlN/C interface, the details at the interface, such as the atomic structure, composition, crystal orientation, and polarity, should be theoretically investigated.
For chemical vapor deposition (CVD) single-crystal diamond, the most important crystallographic surface is the (100) surface, on which most of the homoepitaxial growths and devices fabrication have been performed. Due to the existence of dangling bonds, the (100) surface appears with a 2 × 1-(100) surface reconstruction for the clean surface (C(100)) and 2 × 1-(100):H reconstruction for the surface with hydrogen termination (HC(100)).[11–14] Another feature of H-diamond is the negative electron affinity,[15] which allows the formation of a high interface barrier for diamond heterojunction devices and facile electron emitting source.
In this paper, we systematically investigate the interfacial characteristics of the semi-polar plane AlN(
)/C heterostructures by considering the effects of the terminated planar (N terminated planar (N–AlN), Al terminated planar (Al-AlN)), the surface structure of diamond, and the atomic bonds at the interfaces by using the first-principles study. The different diamond surfaces such as the clean reconstructed 2 × 1-(100) surface (C(100)) and H-terminated diamond 2 × 1-(100) surface (HC(100)) are examined. The effects of the atomic structures of C–N bonding and C–Al bonding at the interface are also analyzed. The density of states, band alignments, charge transfer, and bonding mechanism of the AlN/C are calculated and analyzed. In addition, the N-vacancy (VN) is the most stable native defect in AlN bulk,[16,17] which will bring a change of the carrier concentration when the N is absent in AlN single crystal. Its effect on the band offset at the N–AlN/C(100) heterointerface is also considered in this work.
2. Computational and modeling details2.1. Atomic structure of the interfaceFor diamond, the theoretical equilibrium constant of 3.53 Å was used in the interfacial calculation for consistency. The clean and hydrogen-terminated reconstructed 2 × 1-(100) surface slabs are illustrated in Figs. 1(a) and 1(b), respectively, which are optimized for the calculations. The total density of states (TDOS) of the C(100) reveals the metallic characters, which is mainly related to dangling bond orbital on the dimer atoms, while the TDOS of hydrogenated diamond HC(100) surface still preserves the semiconducting nature. These are in good agreement with the previous computational results.[18,19] For AlN, the theoretical equilibrium value of the lattice constant of 3.08 Å is used. In order to minimize the interface mismatch, the AlN/C hetero-structure slab is obtained by matching an AlN (
) slab to the reconstructed 2 × 1-(100) surface. The relative basis vectors (u′, v′) of AlN(
) surface are extracted through the coordinate transformation of the origin basis vectors (u and v). The process of coordinate transformation is shown in Fig. 1(c). Finally, the slab of AlN(
) surface is shown in Fig. 1(d). Besides, lattice constants x and y at the in-plane correspond to the average of lattice constants for both the AlN(
) slab and the reconstructed diamond 2 × 1-(100) surface, resulting in a 3.0% tensile strain on the diamond side and a 3.0% compressive strain on AlN side.
Four kinds of atomic bonding configurations at the interface are calculated, and configurations are N–AlN(
)/HC(100) (Fig. 2(a)), Al–AlN(
)/HC(100) (Fig. 2(b)), N–AlN(
)/C(100) (Fig. 2(c)), and Al–AlN(
)/(100) (Fig. 2(d)), corresponding to interfaces containing the C–H–N, C–H–Al, C–N, and C–Al constituents, respectively. Here, the surface of diamond is the reconstructed 2 × 1-(100) surface with or without hydrogen, and the AlN(
)/C hetero-structure slab contains 6 layers of C atoms and 4 layers of Al–N.
2.2. MethodsThe total-energy and physical properties were calculated with Vienna abinitio Simulation Package (VASP 5.3.3) based on density-functional theory (DFT) within the generalized gradient approximation (GGA).[20,21] The projector-augmented wave pseudo-potentials and the exchange–correlation functional of Perdew, Burke, and Ernzerhof (PBE) for the GGA were utilized in the simulation.[22] Valence electron configurations for the elements were 2s22p2 for C, 3s23p1 for Al, and 2s22p3 for N. The Monkhorst–Pack scheme was employed for the Brillouin zone integration with a 5 × 9 × 1 k-point mesh.[23] At the same time, the three bottom layers of C atoms and a layer of H atoms were fixed, whereas the others moved freely, in order to facilitate the simulation of a thick diamond substrate and eliminate any strain or stress in the lattice structure as much as possible. After geometry optimization, the interface distance (the distances of C–N, H–N, C–Al, and H–Al along the z-direction) between diamond and AlN was determined by minimizing the total energies of the configurations.[24,25]
The plane wave basis set is included and truncated with a cut-off energy of 500 eV, and the electronic total energy is well converged to 10−6 eV/cell with respect to the k-point sampling. The relaxation of the atomic structure is carried out using the calculated Hellmann–Feynman forces until the force on each relaxed atom is less than 0.01 eV/Å. Finally, a dipole correction is included in the middle of the vacuum to compensate for the spurious electrostatic interactions of the asymmetric slab with a non-vanishing surface-dipole moment.[26–28] Pseudo-hydrogen atoms 3/4H and 5/4H[29] are chosen to passivate the dangling bonds of N and Al atoms at the (
) surface, respectively. Besides, a 20-Å-thick vacuum layer is added to separate the surfaces in adjacent cells in the z-direction to eliminate the interactions between them. As a general rule for calculating the physical properties of the heterostructure, van der Waals correction (vdW-DF)[30] is used to test the effect of the nonbonding forces on interface adhesion formation energy.
3. Results and discussion3.1. Interface adhesion formation energyTo discuss the thermodynamic stability of the four types of heterostructures, the interface adhesion formation energy is calculated according to the following relationship:[31]
where
EAlN/C,
EAlN, and
EC represent the total energies of AlN(
)/C heterostructure, the corresponding part of the AlN(
) and C(100) slab, respectively. The total energies of AlN(
) and C(100) slab are calculated with the same unit cell size, number of
k-points, and cutoff energy as those of AlN(
)/C heterostructures. The calculated interface adhesion formation energies are −1.36, 4.82, −9.21, and −6.34 eV corresponding to the the N–AlN(
)/HC(100), Al–AlN(
)/HC(100), N–AlN(
)/C(100), and Al–AlN(
)/C(100) heterostructures, respectively. After considering the influence of der Waals (vdW) interaction, the corresponding interface adhesion formation energies change into −1.73, 5.12, −9.32, and −6.59 eV. In these heterostructures, the vdW interaction does not alter the interface adhesion formation energy too much. This slight difference between the PBE and vdW-DF results indicates that the vdW interaction plays a minor role in the four examined hetero-interfaces. Therefore, in the following, all electronic properties are calculated without the vdW correction. The positive formation energy indicates that the interface of Al–AlN(
)/HC(100) is not stable, while N–AlN(
)/C(100) is the most stable interface among the four types of heterostructure energetically. Later, we will not show the calculated data for the case of Al–AlN(
)/HC(100) interface. In the present calculations, the order of stability for the four types of AlN/C hetero-structures is as follows: N–AlN(
)/C(100) ˃ Al–AlN(
)/C(100) ˃ N–AlN(
)/HC(100) ˃ Al–AlN(
)/HC(100). The calculation also discloses that the introduction of H at the interface weakens the stability for the AlN(
)/C heterostructure.
3.2. Electronic properties of interfacesThe calculated total densities of states (TDOSs) of the interfaces of N–AlN(
)/HC(100), N–AlN(
)/C(100), and Al–AlN(
)/C(100) heterostructures are shown in Fig. 3. In all these cases, metallic features at the interfaces are observed. Especially, for the Al–AlN(
)/C(100) interface, a dramatic increase of the TDOS near Fermi level is observed. To further understand the spatial evolution of the electronic structure of the AlN(
)/C heterostructure from the interface to the surface, we calculate the projected density of states (PDOS) layer by layer as shown in Figs. 4(a)–4(c) corresponding to the N–AlN(
)/HC(100), N–AlN(
)/C(100), and Al–AlN(
)/C(100) interfaces, respectively. The PDOS is calculated by averaging the total DOS per layer (LTDOS) for each atom in the layer. To clearly identify the electronic states, the PDOSs of the heterostructures are calculated by starting from the interface region and extending to the bulk AlN and the diamond regions. The L denotes the number of the atomic layers and increases with increasing separation from the interface.
Figure 4(a) illustrates the PDOS of the N–AlN(
)/HC(100) heterostructure, in which the interface states mainly result from the interaction between the atoms near the interface. The interface states are contributed from the C 2s2p, and N 2p, Al 3s orbitals as shown Fig. 4(a). It is revealed that the interaction among the atoms at the interface occurs mainly within two C atom layers in diamond for the N–AlN(
)/HC(100) heterostructure, and the electronic properties turn into the feature of the bulk diamond after the third C atom layer. Such a trend in the PDOS varying with increasing atomic layer (L) is similar to those for the other two interfaces of the N–AlN(
)/C(100) and Al–AlN(
)/C(100) interfaces, in which the PDOSs are composed of C 2s2p and N 2p, and C 2s2p, N 2p and Al 3s orbital contribution, respectively. Therefore, for the interfaces of N–AlN(
)/C(100) and Al–AlN(
)/C(100) heterostructures as shown Figs. 4(b) and 4(c), a large number of electronic states appear in the band gap due to the strong interaction between C 2s2p and N 2p or Al 3s near the interface, forming the interface states.
3.3. Charge transferIn order to further disclose the bonding mechanism and the transfer of electrons between the adjacent layers at the interface of the AlN(
)/C heterostructure, plane-averaged electron density difference Δρ(z) along the direction perpendicular to the interface is calculated. The electron density difference Δρ(x,y,z) is described as[32,33]
where
ρ(AlN(
)/C(100)) is the electron density of the AlN(
)/C(100) heterostructure, and
ρ(C(100)) and
ρ(AlN(
)) are the electron densities of the reconstructed diamond 2 × 1-C(100) surface and AlN (
) surface with exactly the same geometries as the AlN(
)/C(100) heterostructure, respectively.
Figures 5(a)–5(c) show the contours of the calculated charge density differences for the interfaces of N–AlN(
)/HC(100), N–AlN(
)/C(100), and Al–AlN(
)/C(100), respectively. Figure 5(a) shows the characteristics of charge transfer for N–AlN(
)/HC(100) heterostructure. It is clearly seen that charge is depleted in the region of Al plane, meanwhile, the charge is accumulated in the region of N plane near the interface. The charge is depleted near the H atoms, which is not shown in the cleaved plane in this calculation. In particular, for the C plane, the charge is depleted near some C atoms and charge is also accumulated near the other C atoms in the first and second layers. Therefore, C–H, C–C, and C–N bonds probably exist at this interface. A similar phenomenon of charge transfer is also observed in N–AlN(
)/C(100) and Al–AlN(
)/C(100) hetero-structures as shown in Figs. 5(b) and 5(c). Comparing the charge transfer in Fig. 5(a) with that in Fig. 5(c), the hydrogen termination on diamond surface weakens the charge transfer from C to N, which is also an important reason for the instability of Al–AlN(
)/HC(100) heterostructure. It is disclosed that the hybridization takes place between C 2s2p and N 2p from the PDOS as shown in Fig. 4(a), and a chemical bonding near the interface is formed, thus stabilizing the interface. The charge transfer only takes place among C, H, N, and Al atoms within three atomic layers near the interface, which indicates a strong localization character.
To further analyze the charge transfer quantitatively, the plane-averaged (x, y) charge density difference Δρ(z) for interfaces of N–AlN(
)/HC(100), N–AlN(
)/C(100), and Al–AlN(
)/C(100) heterostructures are calculated and shown in Figs. 6(a)–6(c), respectively. The magnitude of charge transfer (q) is calculated by integrating Δρ(z) over the full z range. The plane-averaged (x, y) charge density difference for the AlN(
)/C hetero-structure deviates from zero near the interface, suggesting that a considerable number of charge transfer between the diamond and AlN(
) slabs. As shown in Figs. 6(a) and 6(c) for the interfaces of N–AlN(
)/HC(100) and N–AlN(
)/C(100) heterostructures, charge is observed to be depleted noticeably in the C or H and Al atoms the first nearest to the interface, and accumulated in N atoms the first nearest to the interface. In the case of Al–AlN(
)/C(100) hetero-structure as shown in Figs. 6(c), the charge is observed to be depleted noticeably in the C atoms that are the second nearest to the interface Al atoms and accumulated in the C atoms that are the first nearest to the interface. The charge depletion and accumulation are determined by the electronegativity, which is higher for N atoms than for C atoms. It is revealed that the atoms that are the second nearest to the interface contribute to the interfacial bonding. It is observed that the introduction of H plays a role in suppressing the charge transfer between AlN(
) and diamond. The value of charge transfer decreases from 0.037e to 0.025e as shown in Figs. 6(a) and 6(b).
3.4. Band alignmentThe energy band discontinuities or band offsets at the interfaces are also calculated. Usually, there are two methods to evaluate the band offset at the interface. One is the well-established “bulk plus lineup” method[34–36] in which the macroscopic averaged electrostatic potentials are used as reference energies as proposed by van de Walle et al.[37] The other is determined directly from the site projected density of states (PDOS), in which the valence band offset at the interface is the energy difference between the tops of the valence bands of the AlN layer and the diamond layer.[38–40] Although this method can determine the band offset directly, the results are usually sensitive to the computational condition such as the number of k-points, cut-off energy, size of a supercell, etc. Once the alignment of the valence-band edge (EV) and Fermi level (EF) can be obtained, the conduction-band edge (EC) is calculated by the reliable experimental band-gap values.
Here, the band offsets are calculated by the reference macroscopically-averaged electrostatic potential method. In general, the energy of the valence band offset (EVBO) can be expressed as the sum of two terms[34–36]
where the term of Δ
V is related to the lineup of the macroscopically-averaged electrostatic potential, and Δ
EV is the energy difference between the relevant valence band edges (Δ
Ediamond and Δ
VAlN), when the single-particle eigenvalues are measured with respect to macroscopically-averaged electrostatic potential in the buck diamond and the buck AlN, respectively. Due to lattice mismatch between diamond and AlN, Δ
EV is calculated for buck diamond and AlN under a strain of 3% corresponding to the lattice parameters of the hetero-structure slab.
The planar average values of the electrostatic potentials for interfaces of N–AlN(
)/HC(100), N–AlN(
)/C(100), and Al–AlN(
)/C(100) heterostructures are calculated and shown in Figs. 7(a)–7(c), respectively. Each potential is an oscillating function of z that corresponds to the atomic plane. A macroscopic average value of electrostatic potentials as a function of z is calculated as shown in Fig. 7. The calculated results of lineup term (ΔV) and valence band offset (VBO) are displayed in Fig. 7. As shown in Fig. 7(a), the VBO equals 2.0 eV for N–AlN(
)/HC(100) heterostructure, more than doubled VBO value (0.7 eV) of the N–AlN(
)/C(100) heterostructure. However, for the Al–AlN(
)/C(100) heterostructure, the VBO value is 1.5 eV as shown in Fig. 7(c). The terminal surface of AlN has also a nontrivial effect on the valence band offset. The VBO value (0.7) of the N–AlN(
)/C(100) heterostructure with C–N polar interface is less than the VBO value (1.5 eV) of the Al–AlN(
)/C(100) heterostructure with C–Al polar interface. Besides, the hydrogen-termination of diamond also greatly increases the VBO to 2.0 eV from a clean diamond surface VBO value 0.7 eV.
According to the obtained VBOs at the interfaces, the relative position of valence band and conduction-band edges of N–AlN(
)/HC(100), N–AlN(
)/C(100), and Al–AlN(
)/C(100) as obtained from the “bulk plus lineup” method are schematically displayed in Figs. 8(a)–8(c), respectively. Here, the experimental band-gap values of diamond (5.5 eV)[41,42] and AlN (6.2 eV)[43,44] are used to calculate the CBO. The results reveal that the largest band offsets (VBO (2.0 eV), CBO (1.3 eV)) exist at the interface of the N–AlN(
)/HC(100), while the lowest band offsets (VBO (0.7 eV), CBO (0 eV)) are obtained at the interface of the N–AlN(
)/C(100). Each of all the AlN/C hetero-structures shows that it has a typical electronic character of the type-II staggered band configuration. Beside, by gradually reducing the interfacial N content from 100% to 75% at the interface of the N–AlN(
)/C(100), the valence band offsets increase from 0.7 eV for the interface of the N–AlN(
)/C(100) and from 0.9 eV for the interface with VN to 1.0 eV for the interface with 2VN, respectively.
The computed band offsets for AlN(
)/C heterostructures are summarized in Table 1 showing their comparison with the experimental results. According to experimental values of the electron affinity of AlN and diamond,[45,46] the VBO values can also be estimated to range from 0.45 eV to 3.3 eV. In our calculation, the minimum and maximum VBO values are 0.7 and 2.0 eV, respectively, depending on the details of the interface structure. The introduction of hydrogen into the diamond surface dramatically increases the VBO and CBO values. This is in good agreement with the fact that the electron affinity of AlN with Al-termination is higher than that of AlN with N-termination,[47] and also consistent with the negative affinity of the hydrogenated diamond surface. Our data are in discrepancy with the theoretical calculations by Silvestri et al.,[10] where VBO values were 1.6 eV for the interface with C–N constituents and 0.6 eV for the interface with C–Al constituents, respectively. This inconsistency possibly relates to the AlN semipolar surface as well as the selection of the diamond surface structure, which is diamond 2 × 1(100) reconstruction surface in our case and diamond 1 × 1(100) surface in Ref. [10].
Table 1.
Table 1.
![](https://cpb.iphy.ac.cn/html_resources/images/table-icon.gif) Table 1.
Valence band offset (VBO) and conduction band offset (CBO) values calculated by “bulk plus lineup” method. For VBO and CBO, a positive number indicates that valence band top and conduction band minimum of the diamond are higher than those of the AlN.
.
Heterostructure |
Interface atoms |
Formation energy |
VBO |
CBO |
Memo |
Semi-polar AlN(
)/C(100) |
C–H–N |
−1.36 |
2.0 |
1.3 |
|
|
C–H–Al |
4.82 |
– |
– |
not stable |
|
C–N |
−9.21 |
0.7 |
0.0 |
|
|
C–Al |
–6.34 |
1.5 |
0.8 |
|
Ref. [10] |
– |
– |
0.63 to 1.60 |
–0.07 to 0.90 |
theory |
Ref. [2] |
– |
– |
4.0 |
3.5 |
experimental |
Ref. [46] |
– |
– |
2.1 |
1.37 |
experimental |
Ref. [48] |
– |
– |
0.45 |
–0.28 |
experimental |
| Table 1.
Valence band offset (VBO) and conduction band offset (CBO) values calculated by “bulk plus lineup” method. For VBO and CBO, a positive number indicates that valence band top and conduction band minimum of the diamond are higher than those of the AlN.
. |